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We present a new approach to construction of recursion operators for multidimensional integrable sys- 
tems which have a Lax-type representation in terms of a pair of commuting vector fields. It is illustrated 
by the examples of the Manakov-Santini system which is a hyperbolic system in N dependent and N+A 
independent variables, where N is an arbitrary natural number, the six-dimensional generalization of 
the first heavenly equation, the modified heavenly equation, and the dispersionless Hirota equation. 



m ' 1 Introduction 



Existence of an infinite hierarchy of local or nonlocal symmetries usually indicates integrability of a PDE 
system. A standard tool to produce such hierarchies is a recursion operator, see e.g. [3l|29lll2]. In essence, if 
one has an auxiliary linear system and a class of its coefficients for which the inverse scattering problem can 
be solved, then using a recursion operator it is usually possible to write down the most general nonlinear 
system whose solutions can be obtained using the inverse scattering problem in question, cf. e.g. [H 126] . 

Recursion operators can be sought for in a variety of ways, depending on the definition used. For 
instance, a typical recursion operator in (1+1) dimensions is a pseudodifferential operator which maps 
symmetries into symmetries [28] and can be derived from the Lax pair, see e.g. [28l [T7l [LS] and references 
therein. However, these definitions and methods do not immediately extend to higher dimensions [42]. 
Instead, a higher-dimensional equation may admit a bilocal recursion operator as introduced by Fokas and 
Santini, see e.g. [HI [121 [13] , an d [IS]; a prototypical example here is the Kadomtsev-Petviashvili equation. 

From a different perspective, recursion operators are Backlund auto-transformations of a linearized 
equation. They first appeared in this form in several works by Papachristou [301 EH [32] • The idea is that 
symmetries are essentially solutions of the linearized equation and Backlund transformations are the most 
general transformations that relate solutions to solutions. 

However, as pointed out by one of us [23], the same point of view applies to Guthrie's [16] generalized 
recursion operators in (1+1) dimensions. Guthrie introduced them to avoid difficulties connected with the 
lack of rigorous interpretation of action of pseudodifferential operators on symmetries. 

The recursion operators which are auto-Biicklund transformations of linearized systems appear to exist 
only for a certain class of integrable systems, which has yet to be characterized in full generality. Such 
recursion operators are closely related to zero curvature representations whenever the latter exist, see 
e.g. [24J and references therein. 

By a zero curvature representation for a system of PDEs F = we mean a one-form a = ^ Ai dx 1 
which satisfies D x jAi — D x iAj + [Ai, Aj] = on the solution manifold of F = 0; here x l are the independent 
variables and Ai belong to some matrix Lie algebra. Then the operators D x i — adAj commute and we can 
define a matrix pseudopotential \1/ by setting 

V* - [Ai,V] = duAi, (1) 

where du denotes linearization along a symmetry U (see Section 2 below for details). 



For the majority of integrable systems in (1+1) dimensions the pseudopotentials \l/ provide nonlocal 
terms of inverse recursion operators, whereas local terms thereof are, as a rule, limited to zero-order ones 
at most. This approach applies to multidimensional systems whenever a zero curvature representation is 
available. For instance, the recursion operator found by Papachristou in [30] is easily seen to be of this kind. 

Now turn to dispersionless multidimensional systems which can be written as a commutativity condi- 
tion for a pair of first-order linear scalar differential operators with no free terms (i.e., vector fields) and 
no derivatives with respect to the spectral parameter; see e.g. [13J I2TJ HJ [22J [8] and references therein for 
more information on such systems. The systems of this kind are a subject of intense research as they arise 
in a multitude of areas from self-dual gravity, see e.g. [Ml E7J EE1 IE], hyper-Kahler [15] . symplectic [7] and 
conformal [6] geometry to fluid dynamics and related fields, cf. e.g. [351 13 HO]- Even though the recursion 
operators for some of these systems were found, see e.g. [301 ESI I27J [22J [37J [38] , they were obtained using 
either various ad hoc methods or the partner symmetry method, both of which can be applied only under 
fairly restrictive assumptions. 

Below we present a method for finding recursion operators which is based on a generalization of 
pseudopotentials ([1]) using the adjoint representation of the Lie algebra of vector fields. We are convinced 
that our approach applies to a considerably broader class of dispersionless systems than the methods 
mentioned in the previous paragraph. Moreover, our method is also more algorithmic: given a Lax-type 
representation for the system under study, finding a recursion operator of the type described in our paper, 
if it exists, is an essentially algorithmic task while e.g. the partner symmetry approach involves a non- 
algorithmic subproblem of representing the linearized equation as a two-dimensional divergence. 

The paper is organized as follows. In Section 2 we present a general construction of recursion operators 
which are auto-Backlund transformations of linearized systems, and Sections 3 and 4 illustrate its appli- 
cation on the examples of the Manakov-Santini system and the dispersionless Hirota equation. Finally, 
in Section 5 we give a modification of the construction from Section 2 for the case of Hamiltonian vector 
fields and provide some further examples. 

2 The general approach 

Let F = be a system of PDEs in n independent variables x l , i = 1, . . . ,n, for the unknown iV-com- 
ponent vector function u = (w , . . . ,u N ) T , where the superscript 'T" denotes matrix transposition. Denote 
uf in = <9 nH Vln u a /d(x 1 ) 11 ■ ■ ■ d(x n ) ln ; in particular, Wo ... — uOL - As usually in the formal theory of 
PDEs [291 [5j [3], x% an d uf %^ are considered as independent quantities and can be viewed as coordinates 
on an abstract infinite-dimensional space (a jet space). By a local function or a function on the jet space 
we shall mean any function of a finite number of x % and u and their derivatives. We denote 

o N oo „ 

D = — 4- V V v a - 

3 8x i ^ ^—* 'l-'j-lW+^J+l-^^jQi 

a=l ii,...,i n =0 H—in 

the usual total derivatives, which can be naturally viewed as vector fields on the jet space. The condition 
F = along with its differential consequences D l \ ■■■D x £ n F = determines what is called a solution 
manifold, which in general is an infinite-dimensional submanifold of the jet space. In what follows we 
tacitly assume that the total derivatives are restricted on the solution manifold; these restrictions are 
tangent to the latter. 

As usual, the directional derivative along an iV-component vector U = (U 1 , . . . ,U N ) T , where the 
superscript T indicates the transposed matrix, is the vector field 



du a ■ 

a=lil,...,i n =0 ll—ln 



The total derivatives as well as the directional derivative can be applied to (possibly vector or matrix) 
local functions P. Recall [291 E] that U is a (characteristic of a) symmetry for the system F = if U 
satisfies djjF = on the solution manifold. 

Assume now that the system F = can be written as a commutativity condition [X l5 X 2 ] = 0, where 
Xi = Y^j=i£i D x i are vector fields, [• , •] is the usual Lie bracket thereof, and the coefficients £f are local 
functions that may further depend on a spectral parameter A. 

Further consider a vector field Z of the same form, i.e., Z = Yll=i CA^, except that we do not insist 
that (i are local functions. However, we assume that the total derivatives can be extended to £ J , see below. 

The main idea of the present paper is that we look for aniVxn matrix A = (aj) such that 

n 

i=i 
are components of a symmetry U whenever Z satisfies 

[X i ,Z} = d u X l , 2 = 1,2. (2) 

We shall write U = %Ka(U) when such a matrix A = (a") ^ exists; this is precisely a recursion 
operator of the type described in Introduction, i.e., a Backlund auto-transformation for the linearized 
system dyF = 0. Here and below we assume for simplicity that the entries a" of A are local functions; 
however, in principle, nothing prevents them from being nonlocal. 

Note that we do not insist that the vector fields X^ necessarily involve any spectral parameter, but we 
do exclude the case when they involve derivatives with respect to the spectral parameter. 

The condition ([2]) is a system of first-order partial differential equations in the unknowns (i . To show 
that the system is compatible we check the Jacobi identity 

[X x , [X 2 , Z\\ + [X 2 , [Z, X,]} + [Z, [X l7 X 2 ]] = [X h duX 2 ] - [X 2 , duXx) = d u [X 1 ,X 2 ] = 

since [Xx , X 2 ] = is equivalent to F = and U is a symmetry. 

As a rule, the system ([2]) is not solvable in terms of local functions. Therefore, strictly speaking, U 
are not necessarily symmetries of the system F — 0. Instead, they are nonlocal symmetries (or shadows 
in the sense of J2UI IS])- This naturally leads to introduction of pseudopotentials (for instance, ( l and 
their derivatives) and subsequent extension of the total derivatives to include the terms coming from 
pseudopotentials. To simplify notation we shall, however, denote the extended total derivatives by the same 
symbol D x i. Note that when applied to local functions the original and extended total derivatives coincide. 

Note that in a number of examples, where the recursion operators are already known, e.g. the Pavlov 
equation [211 122], our method produces the recursion operators which are inverse to the known ones. More- 
over, the inverses of our recursion operators often have simpler structure of nonlocal terms; in particular, 
this holds for all systems discussed below. Thus, it is often appropriate to invert the operator 9^4 resulting 
from the above construction in order to obtain a simpler recursion operator; the inversion is an algorithmic 
process described in [To] . 

Let us also mention that, in sharp contrast with the case of (l+l)-dimensional systems where one usu- 
ally can make a clear distinction among positive (local) and negative (nonlocal) hierarchies (see, however, 
[2]), the multidimensional hierarchies we have been able to generate contain, an eventual inversion of the 
recursion operator notwithstanding, only a few local symmetries. The same phenomenon occurs for the mul- 
tidimensional hierarchies generated using bilocal recursion operators, see e.g. [23110] and references therein. 



3 The Manakov— Santini system 

Consider the Manakov-Santini system [21] in (N + 4) independent variables x 1 , . . . , x , y 1 , y 2 , z 1 , z 2 and 
N dependent variables u l , 



u y i z2 - u 



N 

; 2z i + Y^ u * u #ti - u i* u Ux>) = o, z = i, . . . , jv. (3) 



As usual, the subscripts refer to partial derivatives. 

System ([3]) can be written [21] as a commutativity condition of the vector fields 



N 



Xt = D yi + \D z i + J2 u k zi D xh , % = 1, 2. (4) 

fe=i 

Assume that Z has the form 

Z = y £ l V j D xi , (5) 

i=i 

(no terms involving D y i and D z k are actually needed). 

It is straightforward to verify that the following assertion holds for any natural N: HU = (U l , . . . , U N ) T 
is a characteristic of symmetry for ([3]) then so is V = (V 1 , . . . , V N ) T , where V^ are determined from the 
equations (|2D with Xj and Z given by (T4j) and (jSJ), that is, 

JV JV 

V$+Al£ + 5>i.l& = X>i^-l&. z = l,...,iV, s = l,2. (6) 

i=i i=i 

To emphasize the dependence on A, the recursion operator given by formula ([6]) will be denoted 9^- 
Applying Di\ to local symmetries yields a highly nonlocal 'negative' hierarchy of the Manakov-Santini 
system. In order to obtain the 'positive' hierarchy with simpler nonlocalities, we look for the inverse ^H^" 1 . 

Inverting the recursion operator *Ha amounts to solving ([6]) for U\. The inverse operator D^ 1 sends 
V = {V 1 , . . . , V N ) T to U = (U 1 , . . . , U N ) T , where U l are determined from the relations 

N N 

Ul a = - V ; s - XV:, -Y,<^ + Y,<^^ i = h--;N, s = l,2. 

3=1 j=l 

Upon multiplying by —1 and removing the trivial contribution XV from U, we end up with the recursion 
operator !lK = — Ji^ 1 — Aid, which no longer depends on A. The components of U = 9t(V) are defined by 
the relations 

N N 

^ = ^+5>^-5>U^ z = l,...,AT, s = l,2. (7) 

3=1 3=1 

The symmetries generated using this recursion operator are complicated nonlocal expressions; the explicit 
form for symmetries obtained by applying ([7|) to the Lie point symmetries of ([3]) is given in the appendix. 

4 Dispersionless Hirota equation 

Consider the equation [HI [6] 

au x u y t + bu y u xt + cu t u xy = 0, a + b + c = 0. (8) 



It has a Lax pair [H] of the form 



\Uy1p X HU t ljj x 

^y = , Wt = , 



u x u_ 

where u = (a + b)\/(a\ + b). 

The vector field Z now can be chosen in the form Z = (D x . An easy computation shows that the 
corresponding recursion operator is given by the formula 

Kx(U) = u x (. 
Here U is a symmetry for flSJ) and ( is defined by the following equations: 

Ct __ Ui_, -f^±\ c UtUx - UxUt ^y.-'hLc ~(^l\ c UyUx - UxUy 



X 



\i u x \u x J x u x A u x \UxJ x u\ 

The inverse recursion operator W = 9 t l^" 1 ([/) is given by the formulas 

TJ/ UyW X UyU X Uy A ~ 1 U X y U t W X U t U X ., , ,, 

Wy = — h — "^ " -U, W t = 1 U. 



UyW X UyU X Uy \~\U X y U t W X U t U X U t fi~lU xt 

y~— + — ^ 7 -U, W t = + 

u x u x A A u x u x u x fi fi u x 



If we replace W by W — U, we obtain a somewhat simpler recursion operator W = 9K(U), where £H 
D^ 1 + Id, so W can be determined from the compatible equations 



Ux A V u x J u x At V U x 



Let us apply *R to the Lie point symmetries, which are X(x)u x , Y(y)u y , T(t)u t , and F(u). To start 
with, we find 91(0) = F(u), where F is an arbitrary smooth function. Upon having agreed to remove this 
trivial contribution from the results, we readily find 

&(Xu x ) = 0, m(Yu y ) = ^^Yu y , &(Tut) = ^—^Tu t . 

A fi 

Thus, Xu x , Yu y , Tu t are eigenvectors of JH. To the best of our knowldege, this is a first known example of 
local eigenvectors for a nontrivial recursion operator. 
Finally, 

X(F(u)) = Fp+^-q, 
where p, q are nonlocal variables defined by the following equations which are compatible by virtue of (|SJ): 

Uy X-lU X y U t U~lU xt Uy A - 1 U t U~ I 

p y = — Px r , Pt = — Px , q y = — q x H : — u y , q t = — q x u t . 

u x X u x u x u u x u x A u x a 

5 The case of Hamiltonian vector fields 

If the Lax pair for the system under study consists of Hamiltonian vector fields, it is natural to apply 
the ideas from Section [2] to the algebra of functions endowed with the Poisson bracket rather than to the 
algebra of vector fields. 

Namely, suppose that the Lax pair for the system under study can (up to the obvious renumbering of 
independent variables) be written as 

Vv- 1 = aip x n-2 + {Hi, i)}, ip xn = bip x n- 3 + {H 2 , tj}}. (9) 



Here {-, ■} denotes the Poisson bracket in question (usually w.r.t. the independent variables x 1 , . . . ,x n ~ A 
only), and a, b are some constants, which are typically proportional to the spectral parameter A. 
Instead of £ J introduce a single nonlocal variable Q defined by the formulas 

n x n-i = aVL x n-2 + {H 1 , fi} + djjHi, VL x n = bVt x n-3 + {H 2 , fi} + dzjH 2 , (10) 

where a, b are some constants which are often proportional to the spectral parameter A. 
Then we shall seek for a recursion operator in the form 

n-2 

U = <R(U) = A Q + J2 A i Q ^ ■ (11) 

i=l 

where now Aj = (ah, . . . , af), j — 0, . . . , n — 2, are iV-component vectors whose entries are local functions. 
If necessary, the terms containing higher-order derivatives can be also included. 

As an example, consider the following six-dimensional generalization of the first heavenly equation, see 
e.g. 0E3]: 

U X p V<i ' ^xt^yz ^xz""yt' \ ) 

It admits [7J a Lax representation of the form (Q, namely, 

1p p = \4>y + {Uy, ifl}, i) q = ~\i) X + {~U X , if}}, 

with the Poisson bracket given by 

{f,g} = D z (f)D t {g)-D t {f)D z (g). 

It is readily verified that (TT2l) possesses a recursion operator of the form U = 9i(U) = Q, where the 
nonlocal variable Q is defined via fllOl) . that is, 

fi p = Xtty + {Uy, Vt] + Uy = XVty + Uy Z tt t ~ Uy t Vt z + Uy, 

Vt q = -X£l x + {—u x , fi} - U x = -X£l x - u xz Vt t + u xt Vt z - U x . 

Upon inversion we obtain a simpler recursion operator U = (D^ -1 + A ld)(U), where U is defined by the 
formulas 

U x = -u xz U t + u xt U z - U q , Uy = -u yz U t + u yt U z + Up. (13) 

For another example, consider the modified heavenly equation [7J, 

U xz = U x yU tt - U xt Uyt, (14) 

which has [7j a Lax representation of the form (151), 

^x = {u x /X, lf>}, 2 = Xipt + {-ut, ip}, 
with the Poisson bracket given by 

{f,g} = D y (f)D t (g)-D t (f)D y (g). 

It is readily seen that (fi~4|) admits a recursion operator U = 9i(U) = Q, where the nonlocal variable Q is 
now defined by the formulas 

Q x = {u x /X, Q} + U x /X = u xy VL t /X - u xt Vly/X + U x /X, 

Vt z = XQ t + {-u t , Q} — U t = XQ t + uu^ly - u yt Vt t - U t . 

6 



Inversion again leads to a simpler recursion operator U = (%K l — A Id)(C/), with U defined by the formulas 

U x = u xt U y - u xy Ut, U t = u u Uy - Uy t U t - U z . (15) 

To the best of our knowledge, the recursion operator ( TT3l) has not yet appeared in the literature, while 
( TT5|) is a special case of the recursion operator for the so-called asymmetric heavenly equation found in [37] 
using the partner symmetry approach. Note that ([TBI also could have been obtained within the partner 
symmetry approach [37]. On the other hand, the recursion operators for the second heavenly and Husain 
equations, which were found in [27J [38] , can be easily recovered using the approach of the present section. 
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Appendix: Symmetries of the Manakov— Santini system 

In this section we use the standard convention on summation over repeated indices. The indices r, s run 
from 1 to 2, the others run from 1 to n. 

The symmetries of the Manakov-Santini system can be routinely computed as solutions U of the 
linearized equation djjF = which must hold only on the solution manifold of (E]); here F l stands for the 
left-hand side of ith equation of ([2]). 

The simplest of symmetries, the Lie point ones, are characterized by the property that U l are linear in 
the first derivatives. A computer-aided computation reveals fourteen Lie point symmetries for ([3]), namely, 

Vli V2, Cl> C2, ", 011, 012, 021, 022, 7, Xl, X2, /£, g<f>, 

where /£, 5 each depend on N arbitrary functions of the coordinates x l , y\ i,j — l,...,N; the left-hand- 
side subscripts indicate the arbitrary functions these symmetries depend on. 
The generators of these Lie point symmetries read 



Vr =u yr, (r=u k z r, a k = y s u k z s, {(5 r s ) k = y 7 'u k s + z r ' u% , j k = u k + y s u k s , 

{x r )k = y r u k + y r y s u k s + ^s^ f f = fkfay^ ^ = -gifayfr*, + ^^ + ^^ Z'. 

This notation is to be read as follows: for instance, the symmetry g <p has a characteristic {g^ 1 , ■ ■ ■ , g 4> N ) 
and the associated evolutionary vector field is J2u=i g^d/dii". 

Let us investigate the action of the recursion operator (jTJ) on the above symmetries. Obviously, equa- 
tions (j7j) determine U k uniquely up to adding the Lie symmetry we denoted f£ (we could also write 
fC, = 5^(0)). Like the integration constants, this term will be omitted in what follows. 

Four classical symmetries, namely £, Ci, C2, «, are mapped to local symmetries again: 

**(/£)=/& H(r)=Vr, JH(a)= 7 - 

The others are mapped to nonlocal symmetries, sharing the same set of nonlocal variables u^ 2 ' k , k = 
1, . . . , N, subject to the equations 



TV 



(2)fe 



Jr +£}«>*'■ (16) 



i=l 



The system f lT6|) is compatible by virtue of the Manakov-Santini system fl3]). Eq. ffT6|) determines a 
covering in the sense of [5], but this covering is infinite-dimensional. Each successive application of £H 
requires one more level of nonlocal variables u^ k , subject to the compatible equations 

N 

u^ k = u^ k + J2<^f. (17) 

Thus we have obtained an infinite hierarchy of successive coverings. 

Upon denoting rp = £H 9 (?aY) = ^ 9+1 \( r ) = Q 9 , we routinely generate 

(l)fc J2)k (2)k k i 

(2)k / .(3)fe (3)fc (2)k j k (2)i . fe j j 

Vr =0 = uy - u K j u l yr - uliuy + U K xi U l xJ U J y r, 

(3)fc JA)k (4)fc (3)fc i <2)k (2)i h (3)i 

T]r = Q = UyT ~ U)J U % y r ~ U xi ' UyT ~ U xi Uyr 

(2)fc j j k (2)i i k i (2)i k i 7 h 

+u x J u x3 u yr + u%u x > Uyr + u k xl u x3 u y J 3 - u k xl u x3 u xh u h yr , 

etc. If we assign level one to the local variables, i.e., u k = u^' k , we observe that the sum in the above 
formulas runs over all homogeneous monomials of the same level. We conjecture that this pattern holds 
for rjr with any natural q. 

In terms of these, we observe the following general formula for (3^ = 9t q (/3): 

{^>f = f^) k + z r Q<f) k . 
Likewise, for jw = 9^(7) = W +1 (a) = a^ q+1 ^ we obtain the expressions 

7 (l)fc = a (2)k = 2M (2)fc _ u i u k^ + yS^ 1 )^ 

7 (2)fc = a (3)k = 3u (3)k _ u m u i _ 2u k iU (2)i _|_ U k xl U xJ ui + y'rfP* , 
7 (3)fc = a (4)k = 4u (4)k _ u m ul _ 2u (2)k u(2)i _ 3^(3)* 

+u x 2 } k u x3 u 3 + u*iuf} l u j + 2u k i u xj u ( - 2 ' )j _ u k xiU l xjW > xhU h + y s -qf )k i 

7 (4)fc = a (5)* = 5M (5)fc _ M (4)^ 4 _ 2M (3)fc u ( 2 )i _ 3M (f M (3)i _ 4u * iU (4)i + u^tfjV? + U^U^V 

+2u x Tu x M 2)j + ul^v? + 2u k xl u x fu^ + Zu k xl u xJ u^ )j ~ u ( x f 'n^u^u' - u k x% u x J V^ 

-u k xl u xj u x g u 9 - 2u k xl u xJ u xg u {2)9 + u k xl u xJ u 3 x gU xh u h + y s r]T . 

Here another pattern can be observed. Namely, the sum again runs over all homogeneous monomials with 
the coefficient at each monomial equal, up to the sign, to the level of the nonlocal variable which is not 
differentiated. It would be interesting to find out whether this pattern holds for all natural q. 
In terms of 7^ we have the following general formula for x = 9^ 9 (x) : 

( x r ) {q)k = y r -f {q)k + z r a (q)k . 

Finally, we have g ^ q) = tH%(j)) = <K (9+1) ( 9 = g ^ q+1 \ 

g ^ k = 9 ^ k = -^ V + |> + u% (-uy + %u> + p*) 



1 / d 2 q k dq k d 2 q k 

+ - ". n . U l v? + 2 -^- u % z s + \ z r z f 

2 \ox l ox J ox % oy s oy r oy s 
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